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Abstract 

The transport of the potential vorticity gradient V<? along surfaces of constant potential temperature 6 
is investigated for the stratified Euler, Navier-Stokes and hydrostatic primitive equations of the oceans and 
atmosphere, in terms of the divergenceless flux vector 13 = VQ(q) x V9, for any smooth function Q of of 
the potential vorticity q. The flux vector B is shown to satisfy a transport equation reminiscent of that for 
magnetic field flux in magnetohydrodynamics. The result may apply to satellite observations of potential 
vorticity and potential temperature at the tropopause. 



14th March 2010 



1 Introduction 



Potential vorticity (PV) is believed to be a particularly significant quantity in the dynamics of the atmosphere 
and the oceans [HE]- For an incompressible fluid the PV density is defined as q = u ■ V# where u> = curl u is 
the vorticity for a fluid with a diverge nee- 1 ess velocity field u and 9 is the potential temperature. This paper 
exposes a mechanism for creating the large gradients in potential vorticity density Vg. Its discussion is based 
on the geometric properties of transport of intersections of level sets of the quantities q and 9 encoded in the 
vector B = VQ(g) x V# which turns out to obey an evolution equation of the forrrj^] 

d t B - curl (U x B) = T> . (1.1) 

Here, the vector field U is formally a transport velocity and will be derived explicitly below in several different 
cases. The quantity T> expresses the rate of change of the flux of B in a frame moving with velocity U. 
When T> = 0, the flux of B is said to be frozen into the flow with velocity ti = u. The Euler equations for 
the stretching and folding of vorticity, and the ideal MHD equations for the frozen-in transport of the flux of 



divergence-free magnetic field, both take the form (1.1) with T> = (see [31 HIE]). However, the vector B 



turns out to be a wise choice even with dissipation : the vorticity dynamics of the Navier-Stokes equations with 



viscosity and the magnetic field evolution in MHD with resistivity, both take the form of equation (1.1) but, 
because T> ^ 0, the corresponding fields are no longer frozen into the flow. Figure 1 shows (model assimilated) 
satellite data for contours of potential temperature 9 on the constant level surface of potential vorticity q = 
2, which lies near the tropopause. See |7] for animations of this data in frames taken every six hours. In these 
animations, the evolution of the contours of 9 on a level set of potential vorticity q = 2 is seen : its appearance 
suggests the stirring of one liquid in another by stretching and folding, such as cream in black coffee. 



Theta on PV=2 




Fi glire 1: A typical snapshot of satellite data taken every six hours of the Northern Hemisphere from ECMWF 7 shows contours of 
potential temperature on a level set of potential vorticity q = 2 located near the tropopause. Notice that high gradients of 9 lie at the sharp 
interfaces of the contours. 



The object of the present work is to derive exact equations of the form (1.1) for the evolution of the 
intersections of level sets of q and 9, in which we will find that T> is given by the divergence-less vector 

T> = - V [qQ'(q) d\vU] x V9 , (1.2) 

for any choice of the smooth function Q, and the vector 14 will be derived in several cases. 



1 The vector 23 = VQ(g) X V# has been discussed in [3] EE] and should not be confused with the cross product V-B X V0 considered 
in [5], in which the scalar B is a steady Bernoulli function. 
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Equations (1.1) and (1.2) are first derived in §2 in the narrower setting of the Euler and Navier-Stokes 



equations. In §3 the corresponding results are derived for the viscous hydrostatic primitive equations (HPE), 
which are commonly used in numerical simulations of the weather, climate and oceans. The stretching and 



folding mechanism inherent in ( 1.1 ) in the context of the viscous HPE may have some application to potentially 



rapid growth of Vq in the atmosphere and oceans, where the occurrence of extreme events is of interest for the 
prediction of variability of the climate. In fact it has recently been shown by Cao and Titi |8] that the solutions 
of the viscous HPE remain regular (see also |9J). It follows that if extreme events do occur in solutions of 
HPE, then these must actually be smooth at sufficiently small scales. Conversely, HPE dynamics, although 
now known to be regular, may still produce extreme events due to the allowed intense stretching and folding 



of and V6> under the dynamics of equations (1.1) and (1.2) 



2 Summary of main results for the Euler and Navier-Stokes equations 

Consider the dimensionless form of the incompressible 3D Euler and Navier-Stokes equations 



Du 



+ 9k = Re~ l Au - Vp. 



D 

TJi 



dt + u-V, 



where the temperature 9(x, t) evolves according to 

D9 



aRe) A9. 



(2.1) 



(2.2) 



Information about V# is needed to determine how 9(x, t) might accumulate into large local concentrations. 
The traditional approach is to study this question through the dynamics of the potential vorticity, defined by 
(u = curl it is the vorticity) 

q:=u-V9. (2.3) 



q = const 




9 = const 



Figure 2 : For the incompressible Euler equations in three dimensions, the vector TS = VQ(g) X V0 is tangent to the curve defined by the 
intersection of the two surfaces q = const and 9 = const. 

The main results for the incompressible Euler and Navier-Stokes equations are summarized in the following: 
Theorem 1 In the cases below, q and 9 satisfy 

d t q + div(qU) =0, d t 6 + U ■ V0 = , (2.4) 
and, with Q(q) as any smooth function of q, the divergence-free flux vector 

B = VQ(q) x V9, (2.5) 

satisfies the stretching relation 



d t B-curl{U xB)=T>. 
The divergence-less vector T> in (2.6) is given by T> = —V(qQ'divU) x V9. 



(2.6) 



1. For the incompressible Euler equations, U = u & thus T> = 0. In this case Dq/Dt = 0, so the 
intersections of the level sets of q and 9 shown in Figure 1 move together with the fluid velocity, u ; 
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2. For the incompressible Navier Stokes equations U is defined as 



u 



-Re' 1 {AuxVe + a^uAO} , q^O. 



(2.7) 



Moreover, for any surface SiTU) moving with the flow ti, one finds 

d 



dt 



dS 



S(u) 



T>-dS. 
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S(u) 



Remark: The three dimensional incompressible Navier-Stokes equations possess only Leray's weak solutions 



while the Euler equations do not even possess these. Thus the manipulations used in deriving (2.4)-(2.8) 
should be considered as purely formal. 



Sketch proof of (2.4)-(2.8) : The derivation of (2.4) follows the same standard manipulations that appear 
in the elegant classic proof of Ertel's Theorem [10J, namely 

Dq (Du _ \ _„ _ (DO 

■ = — oj • Vm • V0 + u) ■ V — 

Dt \Dt J \Dt 

= (Re- 1 Au-V ± 9)-V9 + u-V({aRer 1 A9) 

= d\v(Re- 1 AuxVe + (aRe)- 1 uA6), 



(2.9) 



where V ± 9 = V# x k. The scalar product • V# = and the rest of the terms on the right hand side of 



(2.9) have been regrouped as a divergence. On using divu = one may define U through the equation 

d t q = - 6\v(qu- Re' 1 Au x V0 - {aRe)- l uAe^ =: -div(gW), (2.10) 
in which case divW ^ 0, and 

{dt + U ■ V) 6 = d t 9 + {u- q^Re- 1 [Au xV9 + a^uAO] }-V9 = 0. (2.11) 

The flux J = qti was first introduced by Haynes and Mclntyre in the context of their 'impermeability theorem' 
[HI [12]. There have been objections that U is not a physical velocity [T3 [T4], which have been answered 
by Mclntyre in [TH] but in the context of this paper U has been employed solely as a notational device. The 
remarkably simple form of (2.6) for the incompressible Euler case, in which 14 = u and T> = 0, was derived 



first in [31 El [5] ■ Two versions of the proof of (2.6) are given in the Appendix, the first using Lie derivatives 
and the second using conventional vector identities. ■ 



The right hand side of (2.6) occurs because q is not a scalar function; rather it is a density (a volume 
form). However, because divu = for the incompressible Euler case, it follows that B satisfies 



DB 

~Dt 



BVu 



which is also the standard stretching equation for vorticity u> on replacing 
\B\ 2 satisfies 



1 D 

2D~t 



\B\ 2 = BSB^ A (S) |£| 



(2.12) 

with cj. The squared magnitude 



(2.13) 



where X s (x, t) is an estimate for an eigenvalue of the rate of strain matrix S. Alignment of B with a positive 
(negative) eigenvector of S will produce exponential growth (decay), thus mimicking the stretching mechanism 
that produces the large vorticity intensities that develop locally in turbulence. Ohkitani [16J has studied 
Clebsch-decomposed solutions for u = V/ x Vg where Df /Dt = and Dg/Dt = 0. 

Remark: Moffatt suggested the analogy between the magnetic field in a conducting fluid and the vorticity in 



an incompressible Euler flow [17] (see also [18J). Equation (2.12) continues this analogy. Moffatt's detailed 
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discussion of the topology of magnetic field lines is based on the concept of helicity that requires the existence 
of a vector potential A that satisfies B = curl A where 

a= \[Qve- evQ) + v^. (2.i4) 

The helicity H that results from this definition, 

H= [ A-BdV = [ div (rf)B) dV = <f vpB-ndS, (2.15) 

Jv Jv JdV 

measures the winding number, or knottedness of the lines of the divergence- free vector field B. This helicity 
would vanish for homogeneous boundary conditions. However, if realistic topographies were taken into account 
then the possibility for H ^ would exist. The boundaries may therefore be an important generating source 
for helicity, thus allowing the formation of knots and linkages in the 2?-field. 



3 The case of the hydrostatic primitive equations 

Many simulations of weather, climate and ocean circulation employ the hydrostatic version of the primitive 
equations (denoted HPE). The major difference of HPE from the Navier-Stokes equations lies in the exclusion 
of the vertical velocity component w(x,y,z,t) in the hydrostatic velocity field[^] 

v(x, y, z, t) = (u, v, 0) . (3-1) 

However, this vertical component does appear in the transport velocity field V = (n, v, ew), where e is the 
Rossby number. The velocity field v in (|3.1[) obeys the motion equation 



e(d t + V ■ V)v + k x v + a kQ = eRe^Av- Vp, (3.2) 

and is solved in tandem with the incompressibility condition div V = divu + ew z = 0. The vertical velocity w 
has no evolution equation; it appears only in V • V and is determined from the vertical integral of the incom- 
pressibility condition d'wV = 0. The z-derivative of the pressure field p and the dimensionless temperature 
enter the problem through the hydrostatic equation 

a 6 + p 2 = 0, (3.3) 



which has been been incorporated into (3.2) as its vertical component. The quantity ao is a constant a a = 
H/L 1 is the aspect ratio and R a is the Rayleigh number, which comes from the non-dimensionalization 
of the original equations. By using the vector identity 

V-Vv = -V xC + ±V(u 2 + v 2 ) (3.4) 



in (3.2), the vorticity equation for 

C = curl v (3.5) 

is expressed as 

[d t + V ■ V)C = {aRe^AC + C-VF + curl/ , (3.6) 

where / = — e _1 (fc x v + aofcG). The dimensionless temperature 0, with a specified heat transport term 

h(x, y, z, t), satisfies 

(d t + V ■V)@ = (aRe)- 1 AQ + h. (3.7) 



Equations (3.6) and (3.7) for HPE correspond to those for the Navier-Stokes equations, with the additional 
/i-term. Hence the results in §2 can be lifted over to HPE by defining 

q = C-V6 and B = VQxV6, (3.8) 



2 The primitive equations as used for weather and climate prediction are defined on a corrected spherical grid. 
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where Q(q) can be chosen as any smooth function of the potential vorticity q, which itself obeys the relations 

9 t q + div(qU) =0, q (d t + U • V)6 = . (3.9) 

Here, the formal transport velocity, U, is defined from the vorticity flux density as 

qU = qV - { [J2e _1 AV + /] x V6 + [(ai?e) -1 CAG + /»]}. (3.10) 

Clearly, U includes the effects of rotation within / and the heat transport term h: moreover, divU ^ 0. As in 
the case of the Navier-Stokes equations this is not a physical velocity but is again a convenient mathematical 



device. Accordingly B in (3.8) evolves according to the driven stretching relation 

d t B - curl (U x B) = D (3.11) 

where 

D = -vfqQ'(q)divU) x V0. (3.12) 



This equation is the analogue for HPE of equation (2.4) for the Navier-Stokes case. It implies that the time 



rate of change of the flux of B in (3.8) through any surface S that is transported at formal velocity U is given 
by 

— f B-dS = [ D-dS. (3.13) 
dt 7s(u) Vs(u) 

In particular, when the U-transported surface is chosen to be a level set of temperature (whose normal vector 
is along V0 then the right hand side vanishes, as it should to maintain the tangency of B to such surfaces. 



This means that the effects of PV gradient flux creation due to the right hand side of (3.11) occur only on 
U-transported surfaces that are not temperature iso-surfaces. 

4 Conclusion 

The equations for the evolution of the flux of PV gradient B for Navier-Stokes in §2 and B for HPE in §3 
are the first two main results. Their left hand sides represent the familiar flux transport form that governs the 
stretching processes while the second main feature is the derivation of the right hand side T> and D-terms; 
these diverge nee- 1 ess forcing terms deserve more investigation. Herring, Kerr and Rotunno [19J have performed 



a computational study of vortex re-connection in the Navier-Stokes Boussinesq system (2.1) and (2.2). It is 



possible that the divergence-less vector T> may be the key to understanding this phenomenon with the left hand 



side of (1.1) dominating for early to intermediate times until the effect of T> destroys the frozen-in property. 
A numerical study of the effect of T> may therefore be worthwhile. 

In GFD, topography has been found to have some bearing on the nature of the topology of the B-field, 
because helicity is generated at boundaries, and this potentially leads to the formation of knots and linkages 
in the B-field lines. The atmospheric or oceanic events to which these knots and linkages would correspond 
are not wholly clear. 

Investigation of steady-state balances of B and its interaction with imposed steady coherent shear would 
also be interesting. In fact, analogous investigations of the spatial distribution of the PV flux are already 
underway in other contexts and, for example, have recently been studied via the data analysis of PV fronts at 
the sea surface (Czaja and Hausmann |20J). We hope modern developments in observation and data analysis 
will soon provide new insight into the role and magnitude of the dynamical effects and balance effects caused 
by the transport of PV gradient flux along temperature iso-surfaces. In this regard, see, for instance, the 
recent paper by McWilliams et al [21J for a discussion of filamentary intensification in the ocean by processes 
similar to the stretching of B. Likewise, in the atmosphere, the stretching of B and the associated alignment 
properties of Vg and V# are of interest, particularly in the region of the tropopause, as shown in Figure [JJ 

Finally, the fundamental stretching mechanism in either HPE, or Navier-Stokes is the term B • VU. As 
already noted, the relation B • (B • VU) = B • SB implies that alignment of B along a positive eigenvector of 
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the rate of strain matrix S (noting that divU 7^ 0) will lead to exponential growth in B. Thus, there will be a 
tendency for B to stretch in these positive directions within a large coherent vortex. It has long been observed 
that large-scale vortices develop plateaus in PV and form steep cliff-like edges at the vortex boundary: see 



Rhines and Young [22] and Rhines |23J. As noted above, although the effect of the X>-term in (3.11) is not yet 



clear, it may play a significant role in the mechanism by which plateaus are formed in PV profiles; namely, by 
the transport of PV gradient within a region of nonzero PV and along temperature iso-surfaces. The fact that 
PV gradient flux can be created along temperature iso-surfaces and can penetrate any other U-transported 
surfaces may also help explain the 'leakage', or 'erosion' of PV gradient that is observed in certain regions of 
vortex boundaries. 
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A Derivation of equation (2.6) 



Given the advective transport equations for temperature and potential vorticity, 



DO 

~Dt 



0. 



D 

Dt 



(q d 3 x) = (d t q + u ■ Vq + q div u) d 3 x = , 



(A.l) 



the evolution equation Q2.6| ) for the quantity B = VQ(</) x V0 may be derived easily by using the notation 
of the exterior derivative (d) and the wedge product (A) 



B ■ dS = (VQ(q) x V0) • dS = dQ(q) A d6 . 
The advective time derivative of the leftmost term in this relation yields 

^- (B ■ dS) = [d t B - curl (u x B)] ■ dS along = u 



(A.2) 



(A.3) 



The advective time derivative of the rightmost term in (A.2) yields, using equations (A.l) 

(dQ(q) Ad0 



D 

Dt' 



-d{qQ'd\vu) Ad9 = T>-dS, 



(A.4) 



also along Dx/Dt = u. For incompressible Euler flow, div it = and equation (2.6) arises by equating the 



rightmost terms in (A.3) and (A.4). When u is replaced by 1A for the Navier-Stokes equations, this yields a 



non-zero right hand side because in this case divW 7^ 0. 

The second version of the proof, with the notation ujj = curlW, is simply a direct calculation : 



(VQ) t x (ye) + (vq) x (V0) t 

-V [(qQ' div U) + U ■ VQ)] x (V0) - (VQ) x [V(M • V0)] 

- {V(gQ' divW) + U ■ V(VQ) + (VQ) • VU + (VQ) x u v ) x (V0) 

- (VQ) x {U ■ V(V0) + (V0) • VU + (V0) x w p } 

-V(qQ' divW) x V0 - U ■ VB + (VQ)(cj£/ • V0) - (V9)(uu • VQ) 
+ (V0) x (VQ • VU) - (VQ) x (V0 • VU) 
curl (u x B) — V{qQ' divW) x V0 . 



(A.5) 
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